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Introduction 


Finding  the  roots  of  an  equation  F(x)=0  when  F(x)  involves 
transcendental  functions  and  x  is  complex  usually  involves  some  kind 
of  search  method.  The  efficiency  of  a  search  method  depends  to  a 
certain  extent  on  knowledge  of  the  roots — where  they  are  likely 
to  occur  ir.  the  x  plane,  and  how  many  there  are.  If  a  root  is  knowi 
to  lie  in  a  given  reqicn  in  the  x  plane,  then  a  search  routine  can 
quickly  find  the  root  to  the  desired  accuracy.  But  if  no  information 
about  the  location  of  the  roots  is  available,  a  search  ever  a  wide 
area  must  be  conducted,  and  this  can  be  time  consuming  an!  expensive. 
Consequently,  a  method  for  locating  the  general  area  of  the  roots 
and  determining  the  pattern  of  the  roots  is  very  valuable. 

This  report  describes  a  simple  method  for  graphically  dis¬ 
playing  tne  pattern  of  roots  in  the  complex  plane.  The  advantages 
of  this  method  are  the  same  as  the  advantages  of  computer  graphics 
in  general.  It  provides  a  way  of  visualising  the  system,  in  this 
case,  tile  transcendental  equation-  and  it  allows  one  tc  obtain  a 
qualitative  feel  for  the  Dehavior  cf  the  system  and  therefore 


intuition  about  the  system. 


2 


Determining  the  Root  of  a  Transcendental  Equation 


There  are  a  number  of  ways  to  determine  the  values  of  x  which 
sar.isfy  the  equation  F(x)=0.  (The  case  of  interest  here  is  that  in 
which  F  is  a  complex  function  of  the  complex  variable  x.)  One  way  is 
to  search  for  values  of  x  which  make  Real  (F(x))  and  Imaginary  (F{.:)) 
both  zero  at  the  same  time  where  thc  search  is  based  on  sign  changes 
in  both  quantities.  There  is  a  problem  in  that  a  sign  change  occurs 
across  a  pole  as  well  as  a  zero. 

The  method  used  here  is  that  the  transcendental  equation  is 
written  in  the  form 


f(x)  -  1  (1) 

and  f (x)  is  written  in  the  form 

f(x)  -  R(x)e1?)(x)  (2) 

where  R{x)  is  the  magnitude  of  f{x)  and  4>(x)  is  the  phase.  The  roots 
of  Equation  1  are  values  of  x  which  make 

R(x)  =  1 

*>(x)  -  2mH  m  =  0,  +1,  +2,.., 


Let 


A (x)  -  0.5  if  R (x)  •>_  .1 
h(y.)  0  if  R(x)  <  1 

P  (x)  =  2jH  0  ^  o  <  211 
E(x)  ^  Mr)  +  Rfx) 
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Now  B(x)  is  to  be  the  brightness  on  a  graphics  display.  Since  A(x) 
and  P(x)  eac..  have  discontinuities,  B(x)  will  show  two  kinds  of 
discontinuities.  The  solutions  to  Equation  1  occur  where  the  dis¬ 
continuities  of  A(x)  and  P(x)  occur  simultaneously.  These  points  may 
be  seen  on  a  display  of  B(x).  An  example  is  shown  in  Figure  1 


A  =  0 

A  =  0.5  Jrt  ine  of  discontinuity  ^1  ine  of 
A  tinuity 


/ 


/  , 
/  / 


* 


^  o*J 

2>  \ 
— ' ^  fst.  \  1 


~;>0' 


location  of  a  root 


discon- 
in  <p 


Figure  Is  Illustration  of  a  display  of  B(x). 

A  computer  program  can  be  written  to  compute  and  display  B(x) 
for  a  grid  of  points  in  the  x  plane.  The  roots  of  Equation  1  are 
then  immediately  obvious,  but  more  important,  the  pattern  of  the 
roots  in  the  complex  plane  can  be  seen,  and  a  qualitative  feci  for 
the  transcendental  equation  can  be  obtained. 
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k  Graphics  Display 

The  method  described  above  was  used  to  obtain  the  computer- 
generated  photograph  shown  in  Figure  2.  The  transcendental  equation 


— x 

x^{l+i0.02) 

»  computer  propram  (see  printout  in  Appendix  «  -as  written  to  cal- 
culato  and  display  B(x)  131*131  P—t.  in  the  complex  n  plane 

in  increments  o£  0.1.  Coo'dinate  lines  spaced  3  units  in  x  are 
also  shown.  The  location  of  three  of  the  roots  is  readily  apparent. 

One  of  them,  whrch  is  not  guite  obvious  at  first  glance.  is  almost 
on  the  real  axis  near  1.3.  »»  infinite  n-umbor  of  complex  roots  occur 

upwards  and  downwards  on  the  line  of  discontinuity  in  h.  which  is 
the  line  in  the  general  up  and  down  direction.  The  xines  which  go 

across  the  photograph  are  discontinuities  rn  P.  This  photograph 

.-._a.tr  s  readily  identifiable.  A  search 
shows  how  the  pattern  of  tne  * 

routine  could  be  used  to  find  two  roots  as  accurately  as  desrreu,  ami 
very  efficiently,  because  their  approximate  location  can  be  casLy 

obtained  from  the  photograph. 

The  Photograph  shown  in  Figure  2  was  generated  by  a  computer 
graphics  system  of  the  Computer  Science  nivision  at  the  University  of 
Utah  utilizing  a  multiple  point-by-point  exposure.  »  Vagram  of 
the  system  is  shown  in  Figure  3.  The  approximate  Univae  1100  empu-er 
time  used  to  make  the  photograph  was  HO  seconds. 


* 
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x 


r 


3  units  in  x/di vision 


Figure  2:  A  computer-generated  photograph  showing  B(x)  in  the  complex 
x  plane.  There  are  131x131  points  corresponding  to  increments  in  x  of 
0.1.  The  roots  occur  at  the  intersections  of  the  lines  of  discontinuity. 
Tnree  roots  are  shown  in  this  photograph.  Ore  of  them  is  almost  on  the 


real  axis. 
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Poor  Man’s  Graphic  Mot hod 


All  th.  in  formation  contained  in  Figure  2  in  also  contained 

in  the  printout  shown  in  Figure  3,  even  though  the  printout  is  not 

as  pretty.  The  direction  or  real  x  i  <=  rf,„„vc  j  c  i 

j.ear  x  js  reversed  from  that  of  the 

Photooranh.  The  safe  basic  fcttod  Ms  uscd  to  ob,,in  thc  pci„tout 
m  Figur..  *,  but  the  program  is  much  rlmpi  *r  because  no  graphics 
system  is  involved,  and  the  method  can  be  used  with  any  digital  com- 
P""°3  "yStem*  Tnc  procT‘ait;  listing  is  given  in  Appendix  B. 

In  the  printout  shown  in  Figure  3  each  space  represents  a 

^  pl~X  /alU  51  "  ai'd  tne  lncreraeniS  i* led  here  are  0.1.  7.  slash 

in  a  space  indicates  p (  >- >  i  ,  .  •  .  . 

-  1'  ifle  value  of  6(x)  are  in¬ 

dicated  «S  : 


Blank 

Period 


Xn  “‘1R  Valur  of  *  was  calculated  in  floating  point  and  1  ion 
convort'prt  i-rt  M  v/  A 

- -v  *.!,».  i  point,  a  2  moans  2.0  <  a  *t  ■>  »  i 

_  J-u-  Thj  locus  of  points 

where  $  chances  from  2P  to  o  * rri'ii 

a‘~  carxly  apparent,  and  the  other  gym- 

bols  show  hoc  4  is  changing  over  the  piano,  since  the  prog™  is 
crtro  ely  si,,, pie,  this  oethcd  should  hu  ,  very  useful  tool,  even  to 


tiioso  without  access  t 


nics  .system.  The  I1O8  computer  time 


U  °d  t'1  0,'tair‘  the  print°uu  f  Pi9ur«  3  was  13  seconds.  Tne  method 
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Ficure  3:  A  printout  showing  the  san»e  information  ar  the  photograph 
in  Tiquro 
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was  used  hero  at  the  University  of  Utah  with  good  success  in  solving 
uhe  dispersion  equation  describing  the  double-stream  microwave 
interaction  between  lectrons  and  holes  A  semiconductors . ’  A 
litfe  experimentation  with  density  of  characters  would  result 
in  a  more  pleasing  printout. 


Conclusions 


The  inherent  power  of  computer  graphics  has  been  applied  to 
the  solving  of  transcendental  equations.  Th^.  method  described  here 
is  simple,  yet  powerful.  It  provides  a  way  to  v:sualizo  the  character¬ 
istics  of  the  transcendental  equation  and  obtain  insight  into  tne 
nature  of  the  roots.  W.th  this  insight,  a  search  routine  could  be 
used  very  efficiently  to  obtain  the  roots  to  the  desired  accuracy. 


5 Lewis  C.  Goodrich.  A  small -signal  field  analysis  of  double¬ 
stream  interactions  in  finite  semiconductors.  Ph.D.  thesis. 


University  of  Utah,  Salt  Lake  City,  Utah,  June  1970. 
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APPEND  IT!  A 


This  appendix  contains  the  Listing  of  the  program  used  to 
obtain 'the  photograph  in  Figure  2. 
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APPENDIX  B 


This  appendix  contains  tile  listing  of  the  program  used  to 
obtain  the  printout  given  in  figure  3. 


001  dl 

i  *- 

Ch,  put  X » t 

0Cii,3 

p  •*. 

Olf  tfiSIOU  ir  ( 1 3 J  )  ,  I A  ( 1 

OCll.4 

» 

DAT  A  K/lH  »  1H,»  1  HP » 

bf)  j  1.0 

4  «• 

Da  1  A  /l!i//  »  M3/1M 

uGl 1 1 

Cm. 4.  ij  liKlt’ 

UOi  u 

u  4 

Xh-- 7 . 

00113 

n 

VI-  f. 

00114 

1  f- 

nxi--. i 

001  ib 

nx*<=  •  I 

001  lt> 

DO  u  u“i»i31 

UUlPl 

1  1  + 

Du  1  “J. >1.1 

UUi,  0 

i<-  » 

X~t.Kr'i_  l  Xh  »  /.  L  ) 

001  P~> 

In*. 

F-i  ./  (tcXJ  (/.)  ♦  <y-<  1  .  »  . 

0  0  i  p  b 

lM 

A:  ,  ( P  } 

0  fl  1  r.  / 

lb* 

Ji-  (  t .  1  .  J  1  a  ( I  )  -  , 

GO  131 

1..‘ 

IF  (n  -1.L1  .J  .)  iAU)-'  3 

00130 

1  /  * 

P  —  ( r  ) 

00134 

lor- 

A I  —  /tli  it,  ( f- ) 

0  0 1 0  0 

I  M  1 

I  *  i  i  J  —  A 1  • .  i  j  2  ( A I  » K  ) 

HJl.'.O 

<_l;  t 

It  (ht.l  .LI  .il.  )i-»HI-P>»T  t 

UUI  40 

P  i  ; 

K ~Fi  ii 

001-41 

CP  r 

Id ( i )  -  r c  (;% + 1 ) 

UOJrt*: 

P  'i  > 

i  i-  dx* 

UO  i  4  4 

P  ‘4  t 

DKIi.T  flP 

l  0  ipt: 

P‘,' 

3  FOi.ihAl  U/.  131  a  1  ) 

nt)  lbO 

fcf  >  1 

Pit  Ifil 

UOii.i 

t-  /  * 

n  Fi*kj*«a  l  (  x H  i  t  1  A 1  ) 

00 1  nt 

Or. 

Xi.=-7. 

001.0 

PO-t 

a  y  i-xn  x i 

t  *  U 1 !  L> 

ji,  i 

STt.D 

unit  o 

31  * 

I”  1 

>1C  7) 

*  t  H'J » 


i))  ) 


lh- 


1H 


